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We consider the pairing between conduction band electrons, and the valence band holes in the 
neutral bilayer-type structures. By employing the bilayer Hubbard model, we show the possibility 
of the inter-plane exciton formation in the system without applied external field. The in-plane and 
inter-plane Coulomb interaction effects on the pairing mechanism are considered, and the role of 
the in-plane particle hopping asymmetry on the gap behavior is analyzed in the paper. We show 
that both Frenkel-type pairing channel and Wannier-Mott-type excitonic pairings are present in the 
considered system. We analyze also the structure of the chemical potential in the bilayer system. 
The temperature effects, and the tunable inter-plane electron hopping effects are discussed. For 
the Frenkel channel, we have shown a particular behavior of the chemical potential at the very low 
temperatures, which is related to the degenerated Frenkel-gap. 


I. INTRODUCTION 

The physics related to bilayer structures is now the 
subject of intensive theoretical and experimental re¬ 
searches due to the recent developments in the field 
of electronic nanotechnologies and researches on the 
high speed information processing technologies. The 
most prominent example in the technological applications 
of such systems is related to the recently investigated 
graphene-based bilayers, which are very promising in the 
context of future nanoelectronic devices i - — Another in¬ 
terest related to the double layer structured materials is 
that they are good candidates for the exciton conden¬ 
sate formation at the very low temperatures*^ - — The 
possibility of formation of the excitonic condensates in 
semiconducting materials has been suggested long years 
before22r— and, it was an intensively studied problem in 
the intervening years 

To control the formation of the excitonic condensate 
in the double layer structured systems, an external elec¬ 
tric or magnetic field are supposed to be applied to the 
bilayer structure. The possibility of the excitonic con¬ 
densate formation in two-dimensions (2D) , under strong 
magnetic field in a semiconductor heterostructure, with 
a very field-sensitive band structure, is given in Refs. |32| 
and ^3. A thermodynamically stable condensate forma¬ 
tion in InAs-GaSb -based system is discussed in Refill, 
where a study of the Keldysh-BCS theory is given. Mean¬ 
while, a new platform-material is proposed for studying 
the excitonic condensates. Recently, the detailed phase 
diagram for the inverted InAs/GaSb bilayer quantum 
wells is given also^ and the possibility of the stable s- 
wave exciton condensed phase is proved. Contrary, for 
the large interlayer tunneling amplitude, a topologically 
non-trivial quantum spin Hall insulator phase is obtained 
with the p-wave pairing gap. 

Despite of many theoretical investigations in this field 
of research the exact results, for the ground state prop¬ 
erties of the electron-hole systems are rare. A determi¬ 
nant quantum Monte Carlo study of the exciton conden¬ 
sation in the bilayer Hubbard model is given recently*^ 


Both, on-site Coulomb interaction and interlayer electro¬ 
static interactions are considered.However, the question 
whether the true excitonic condensate formation is pos¬ 
sible at the very low temperatures is still open. 

Another interesting problem in the field of the cor¬ 
related bilayer systems^ is related to the phase coher¬ 
ence in the excitonic condensate regime. Mainly, it has 
been shown recently that the excitonic pair formation 
and the excitonic condensation, are two different phases 
of matter j 30 i 31 i 35 i 36 Thus, the question of the coherent 
excitonic condensate formation in double layer structures 
needs to be revised substantially. 

In this paper we examine the possibility of the exciton 
formation in the asymmetric double bilayer electron-hole 
structure with zero applied external electric field. We 
suppose the case of the negative hopping amplitude for 
the layer-1 electrons (see in Fig. [Tj) , rending the problem 
to that of the electron-hole type bilayer structure. We 
will show how the excitonic gap is formed and how the 
strength of it determines the type of the excitons, which 
are divided into two categories: Frenkel and Wannier- 
Mott-type excitons. Frenkel-type excitons are mediated 
by the localized Coulomb interaction and are the strongly 
bound excitations, and the Wannier-Mott excitons are 
delocalized, controlled by a weak Coulomb attraction be¬ 
tween the electron-hole particles. For a asymmetric type 
bilayer structure we demonstrate that the solutions of the 
chemical potential appear in the form of two well defined 
bands, separated by a gap region and we show that the 
values of energies, at which the Frenkel-type formations 
could appear in the system, are located between differ¬ 
ent bounds of the chemical potential forming the Frenkel 
channel in the exciton formation energy scheme. In the 
region of the very low pairing gap values, we show the 
presence of the Wannier-Mott-type excitonic gap lobes, 
which decrease when increasing the intralayer Coulomb 
interaction. 

We show, that at zero temperature the excitonic 
gap parameter is degenerated, which is a direct conse¬ 
quence of the presence of the Frenkel-type energy zone 
in the chemical potential solution and which disappears 
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at higher temperatures, when the single valued Frenkel 
channels are opening. We will calculate the temperature 
dependence of the excitonic gap parameter by consider¬ 
ing different values of the intralayer Coulomb interaction 
parameter U/\t±\. 

The paper is organized as follows: In the Section [TT| 
we introduce the bilayer Hamiltonian to describe the 
possible interactions in the system. We discuss also 
the linearization of the non-linear interaction terms in 
the Hamiltonian and the general solution strategy of 
the problem. In the Section [III] we introduce the effec¬ 
tive fermionic action and we introduce the self-consistent 
(SC) equations for the excitonic gap parameter and 
chemical potential. In Section El we give the numerical 
solution of the problem of the interacting bilayer system, 
and we present the results for the excitonic gap param¬ 
eter and the chemical potential. The temperature de¬ 
pendences of calculated quantities are also given there. 
Finally, in the Section [V] we give the conclusion to ob¬ 
tained results. 


II. THE METHOD 

We introduce the Hamiltonian of our bilayer system 
with the in-plane and inter-plane couplings and we dis¬ 
cuss the theoretical formalism. The bilayer is composed 
of two square lattices numbered as / = 1,2, doped respec¬ 
tively with electrons and holes (as it is shown in Fig. [Tj) . 
The Hamiltonian is 

*=-£*« ( C L C J»* + h - C •) (4l<T C i2a + h.C.) 

(i,j) i,c T 

l,a 

- ^2 + u ^2 [fa»/f - 1/2) • (rii 4 - 1/2) - 1/4] 

i,l,cr i,l 

+W J2 [("«* - 1/2) • (n i2rj , - 1/2) - 1/4]. (1) 

i,<j,cr' 

Here, (i,j) denotes the sum over the nearest neighbors 
lattice sites, ti,l = 1,2 is the in-plane hopping ampli¬ 
tude and t _l is the inter-plane hopping amplitude. For a 
simple treatment, we consider the double layer structure 
without applied external electric field and we suppose 
that the chemical potentials in both layers are opposite 
in sign [i i = —/i 2 - This, in turn, will introduce the p- 
doping in one layer and n-doping into another. The on¬ 
site Hubbard -U term in Eq. (jTj) is written in the way that 
the case fi — 0 corresponds to the half-filling. The pa¬ 
rameter U is the local intralayer Coulomb interaction. 
We consider the double layer with respect to the half¬ 
filling conditions in each layer (rq) = 1, for l = 1,2. The 
formation of the excitons in the layered structure and the 
possibility of their further condensation at the low tem¬ 
peratures requires the attraction between the electrons 
and holes. In order to describe the interaction between 
the layers in the structure, we include also the Hubbard 


n-doped layer-2 
1 - th lattice site 



p-doped layer-1 

FIG. 1: (Color online) The electron-hole bilayer structure, p- 
doped hole layer (see the layer-1) and n-doped electron layer 
(see the layer-2) are represented. The outermost diagonal 
sites in the layer-1 represent two-processes: the electron hop¬ 
ping, from layer-1 to layer-2 (see at the left diagonal lattice 
site) and the electron-electron interlayer interaction, before 
hopping process (see at the right diagonal lattice site). 


interaction term between different layers and this is de¬ 
scribed by the last term in Eq. ©• The parameter IF, in 
Eq. (jTj) , is the local interlayer electron-electron Coulomb 
repulsion. This term will be responsible for the excitonic 
band formation. 


A. Decoupling of interactions 

In order to decouple the quadratic terms appearing in 
the expression of the bilayer Hubbard Hamiltonian we 
will write the terms in the Eq. m in more convenient 
form, by shifting the chemical potential in both layers 
and by adding to it the interaction constants. Namely, 
we introduce /q, Z = 1,2, as /q = tn + U/ 2 + W. The in¬ 
tralayer Coulomb interaction term could be also written 
in more symmetric form, taking into account the spin- 
rotational invariance, namely 

71 ? 

n in n U i = - Sf z (2) 

with mi = rni\ + nui and 

Siz = - (jiu^ - nui ). ( 3 ) 

Next, we will pass to the Grassmann representation 
for the fermionic variables and we write the partition 
function of the system by employing the imaginary-time 
fermionic path-integral method. 3 - For this, we introduce 
the imaginary-time variables r at each lattice site i. The 
time variables r vary in the interval (0, /J), where (3 — 1/T 
with T being the temperature. The time-dependent 
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variables cu(t) (c^(r)) are satisfying the anti-periodic 
boundary conditions for fermions cu(r + /3) = — cu(r)m — 
Then, the grand canonical partition function of the sys¬ 
tem is 


Z = 



o-s[^ A 


(4) 


where, the action in the exponential is expressed as 

rd 

S [ct,c] = Sb [c+,c] + / drH (r). (5) 

J o 

The first term, in Eq.(j5j), is the fermionic Berry-term. It 
is given as 


Sb 



Furthermore, we will combine the quadratic and linear 
total density terms in Eq. © and we will decouple the 
obtained non-linear term by the scalar-field Hubbard- 
Stratanovich linearization procedure^ 7 For the total elec¬ 
tron density quadratic term, we have 


exp 



= J [DV] exp 


y_2 

E — ij +iVa i nu 



(7) 


next, for the electron density difference quadratic term, 
we have 


exp 


U_ 

4 


E ( n ^t - n m ) 2 

i,l 



A lil 

u 


+ A ci / • (nu^ - riui) 


( 8 ) 


We have introduced in Eqs. © and © two real field vari¬ 
ables V and A c , which couple to the total density and 
density difference terms, figuring in the initial Hamilto¬ 
nian. After those transformations, the partition function 
of the system takes the following form: 

Z = j [DYDc] J [DV] J [DA c ]e' s# l ct ’ c ’ y ' 4 'l, (9) 


where, the action in the exponential is the total fermionic 


action of our system. It is given by 
So [c f ,c, V, A c ] = Sb [c f ,c] 

+ Yl dl~ti {c\la( T )Cjla(r) + h.c?j 


(hi)' 

l,a 


st± u dT ( 4 i <7 ( r ) c * 2 <T (r) +h.c. 




E / dr 


0 r K 2 i(r) 


U 


-iVu{T) • ( nu(T ) - 


u 


0 TA2 

dr 


i,l 


-W 


A^(r) 

U 


- A c *i(r) • (n i/t (r) - n i 4 (r)) 


rf' 

Li 


dTn ilcr (T) ■ n i 2 a '(T). 


( 10 ) 


For the integration over the Vu -variables, we have the 
following scalar field integral 


/[W] 


D -E;,z/(f dr\S^-iV iV {nii{T )-^)j 


= const • e 




( 11 ) 


where the action in the exponential in the right hand side 
in Ea. (pT]) is given after the saddle point evaluations for 
Vu. It is not difficult to see that 


Vi 


s 


i_U_ 

T 



( 12 ) 


with fii = (rq), l = 1,2, being the quantum statistical 
fermion averages for the total particle numbers in each 
layer. Then the action S [F 5 ] in Ea. ([TT|) has the form 



The integration over the scalar field A c u could be also 
done by the saddle point method. For the saddle point 
value of A ci i we get 

Af* = | <S«*(r)> (14) 

and the corresponding action becomes 

s [Af] = - E / dT {- l Y^ S d n iicr- (15) 

• 7 „ 7o 


For the half-filling case, we will fix Af x = Af 2 = 17/2. For 
further simplifications, we will write the interlayer inter¬ 
action term in Eq. m in a compact form, by introducing 
the following complex variables 

A i,aa'(r) = 4 ia( T ) C l^'( T ) 


(16) 
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with all possible spin directions and their complex con¬ 
jugates. The variables A^ acr /(r) are linear in electron 
density. Then, the interlayer Coulomb interaction term 
will reads as 


S [At, A] = -W V f dr\^{T)\ 2 . (17) 

/ J 0 


This form of the interaction action is more convenient for 
further decoupling process. 


where the action contains already the terms given in 
Ea.(H9l): 


S [c^c,Q\Q] = S B [c\c] + 

+ F f dTti (4 + h - c •) 

/. .v Jo v 7 

(hj) 

i,<j 

r 0 

+t± W / dr 

J o 


( c lla( T ) C *2a('T) + h.C. 


+S [V s ] + 5 [Af] + 5 [At, A, Qt, Q] . 


( 21 ) 


B. The interlayer interaction 


In the next Section, we will use the form of the action 
given in Eq. m to derive the expression for the excitonic 
gap parameter. 


We can decouple the non-linear term in Eq. (H7l) by in¬ 
troducing the complex Hubbard-Stratanovich decoupling 
fields Q\ gg i (t) and coupled to the complex lin¬ 

ear variables Ai j<7<7 /(r) and AJ acr , (t). The transforma¬ 
tion is 


exp 


W 


E / dr \ Ai ’° 

i rr rr I " 0 


J [DQ^DQ]exp - ^ J dT^\Q it<7<7 >(T)\ 


t 


+Qi,aa' (t) + A]^, ( T)Qi^ a , (t) 


(18) 


III. EFFECTIVE ACTION FOR FERMIONS 


Now, it is easy to write the total action of the system 
given in Eo. ([ 2 T|) in the Fourier transformed form, using 
the expression of the decoupled action in the Eo. ([T9j) . 
To this end, we will introduce here the following four 
component Nambu spinor Tk (v n ) 


K) = 


^ki -\{yn) 

C ki4,(^) 
c k 2 
c + 

C k 2 t 


On) _ 


( 22 ) 


and the complex conjugate Nambu vector Here 

k, is the electron wave vector in the reciprocal Fourier 
space, and the fermionic Matsubara frequencies^ v n = 
(2 n + 1 ) 7 r//3, with n = 0, ±1, ±2,.... Then, we have 


Furthermore, we will see that the decoupling fields 
Qiaa'( r ) an d Qi,cra' (j~ ) play the role of excitonic pair¬ 
ing gap parameters for different spin orientations. The 
action related to the interlayer interaction term will reads 
as 


S 


[At,A,Qt,Q] = W 

/ JO 

i rr rr' 


P dT \\Qi,aAr )\ 2 


W 


+QLa> ( T ) A*,™' (r) + Aj , (r)Q i)<T<T , (r) 1 . (19) 


Thus, the last term in Eq. (p~Q|) will transforms as 
Sw [ c \ c ] ~^ S [A - * - , A, Qt, Qj. Next, we will write the 
partition function of the system 


Z = J [Dc'Dc] J [DQ ] DQ] e -S[c\c,Q\Q}^ 

( 20 ) 


Z = J [D^DV] J [DQ^DQ] e -5[» t ,»,Q t ,Q] ) 


(23) 

where the action S [4^, \Eq Q] in the exponential is 
5 [&, ¥, Q\ Q] = ~^J2 *!> n)G? K) $kK). 

k ,V n 

(24) 

Here, the matrix G^ 1 (v n ) is a 4 x 4 matrix, given as 

K) = 

( Elk {yn) 0 ^_L T Qf\- 0 ^ 

_ 0 -Elk (. Vn ) 0 -t_L - Q{ ; 

t± + Q{ t o F 2k {yn) o 

\ 0 —t _L — Q 44 , 0 —F 2 k {yn) / 

(25) 
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The energy parameters £)k (^ n ), l = 1,2, entering into 
the Eg. ([25]) , are given by 

Ei\a (v n ) = —iis n — Af + + fli - 2 ^". (26) 

Next, t/k, Z=l,2, are the renormalized in-plane hopping 
amplitudes ti k = 2 £/ 7 k, and the energy dispersion is 

7 k = cos(k x ) + cos (fey). (27) 

It is not difficult to derive the saddle point expression for 
the decoupling fields Qi, a a'- Indeed, from Eqs. (fl9|) and 
{UD, we find that 

Qaa' = -W (c\ ia (r) C lia , (r)) , (28) 

and 

Qi a , = -w[c\ ia (r)c 2i<7/ (r)). (29) 


We have 

Qua’ = Ql a t = -W (Ko') 

= -JBNY E ( C 2kdr(* y «) C lk<r , (* y n)) • (32) 

^ ' k,zy n 

We will rename here the gap parameters, as Q aa > = A acr /. 
From the structure of the inverse Green function matrix 
given in Eq.(j25j), it follows that A acr / = A aa 5 acr r . We 
should mention here that this result is true only for the 
bilayer without applied external electric field, such as in 
our consideration. For the general case, when an external 
electric field is present, then A|* 7 ^ AJ* and \fi\\ 7 ^ |/^ 21 - 
We will not consider here this case and it will be the 
subject for our future considerations. After calculating 
average in Eq. ([32]) and taking into account the half-filling 
condition for the total particle densities in each layer, we 
have found the following system of coupled SC, non-linear 
equations 


The quantum statistical averages in Eas. ([28|) and ([29]) 
are given with the help of the fermionic action in Eq. , 
mainly 

(...) = Z~ x J [£W t .D’P] J [DQ ] DQ] 

(30) 


1 

N 


J2 np ( e k) + n F (e k ) 

k 


= 1 , 


(33) 



(£_i_ + A) • [rip ( 4 ) ( e k)] 

y/F* + (t ± +A ) 2 


(34) 


In fact, the parameter Q aa ' in Eq. ([28]) plays the role of 
the excitonic gap, which is formed between the p-doped 
hole layer-1 and n-doped electron layer-2 (see in Fig. [Tj) . 
Particularly, the gap Qtt describes the pairing gap be¬ 
tween the electron with the spin j" in the layer -2 
and the cikt = d\^ w ^h the spin j. in the layer -1 
(we denoted by d \the hole creation Grassmann opera¬ 
tor). This corresponds to the conventional excitonic pair¬ 
ing picture as in the usual bulk semiconductors j ^ ~ 29 ! 35 ! 36 
Contrary, the gap describes the pairing between the 
electron with spin f and the cik^ = d\^ hole with 
the spin This is the result of the excitonic pairing with 
the possibility of the spin-reverse, when the electron is 
created in the layer-2. This type of pairing mechanism 
is rare, because of the fast e-h recombination effects that 
are always present in the considered system. For simplic¬ 
ity, we will not consider this type of pairing mechanism 
for the first treatment. 

Without any loss of generality, we suppose the 
case of the pairing states with the uniform real gap 
parameters ] 35 ! 36 Remember, that the interlayer pairing 
fields Q t and Q are local in space and the part of the 
action corresponding them is then 

Sp [Q\Q,c\c]=-^Y,Y,^’ x 

kjZ'n CT,<j' 

l,V 

X (1 - ^/') c Lo- C k/V'- 

( 31 ) 


where rip (z) is the Fermi-Dirac distribution function 
rip (z) = \/[eP z + 1 ) and the parameter Fk is 

tik ~ ? 2 k - F (35) 

For the convenience, we have omitted the spin indexes 
near the excitonic gap parameter A. Furthermore, the 
parameters e k , are given by the following relations 

^k = — "7 + 2 (t\ + £ 2 ) £k + W ± F k + (t _l + A) 2 . 

(36) 

The chemical potential /i is determined so as to maintain 
the average particle density at half-filling and this is rep¬ 
resented by the first equation Eq. ([33]) in the system of 
SC equations. 


IV. THE RESULTS 
A. Gap solution 

In this Section we present the numerical solution of 
the system of coupled SC equations given in Eqs. ([33]) 
and We will perform the solution of this sys¬ 

tem of equations by employing the finite-difference ap¬ 
proximation method in which the fast-convergent New¬ 
ton’s method^ is employed for the nonlinear equations. 
The accuracy of convergence for numerical solutions is 
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achieved with a relative error of order of 10 -7 . The k- 
summations in Eas. ([33|) and (IM1) were performed with 
the (100 x 100) ^-points in the First Brillouin Zone (FBZ) 
on a simple square lattice layers. In Fig. [2] we have pre¬ 
sented the band energy dispersion curves for the inter¬ 
acting bilayer system. The calculation was performed 
for the fixed value of the intralayer Coulomb interaction 
U/\ti \ = 1.6. The interlayer Coulomb interaction is cho¬ 
sen W/|ti| = 0.75 with the interlayer hopping amplitude 
t_ l = 0.04|£i| and we have sampled the FBZ by a number 
of 100 ^-points. The functions in Eq. ([36l) are represented 
on a discrete equally-spaced mesh of the sample square 
lattice. We see in In Fig. [2j that very small charge gaps 
AG = eik — C 2 k are opening in the system along the high 
symmetry directions: (0,0) —)> ( 0 , 7 r) (at the boundary, 
the region numbered by 1), ( 0 , 7 r) —> ( 7 r, 7 r) (see the re¬ 
gion numbered by 2 ) and (tt^tt) ( 0 , 0 ) (see the region 
numbered by 3). The inset, in Fig. [2j shows the high 
symmetry points on a 2D square lattice. In Fig. [3] we 
have shown the Monkhorst-Pack grids for band struc¬ 
ture presented in Fig. [2j The left panels in the picture 
represent the Monkhorst-Pack scheme, when the band 
structure grid-points are connected. The panels from top 
to bottom represent the Monkhorst-Pack sets when eval¬ 
uating from: 7 —>> X (the left-top a panel in Fig. [2j), 
X M (the left-middle panel in Fig. [2]) and M —>> T 
(the left-bottom panel in Fig. [2}. The right panels in 
Fig. [2] show the Monkhorst-Pack sets, when the grid- 
points are unconnected. The same high symmetry di¬ 
rections are considered when evaluating from top-right to 
bottom-right panels. We see clearly in all panels in Fig. [ 2 ] 
(especially it is visible in the panels with the connected 
grid-points) that a charge-gap like small gap is always 
present in the band structure and it is diminishing, when 
evaluating along high symmetry directions clockwise (i.e. 

r^x^M^r). 

In Fig. |4] we have presented the [/-dependence of the 
excitonic gap parameter A for different values of the 
hopping-normalized temperature T/\ti\. For the energy 
parameters in Fig. 0J we have considered the follow¬ 
ing values W/|£i| = 0.4, t\ = —1.0, 1 2 = 0.5|£i| and 
t± = 0.021ti |. The |ti | energy-parameter serves as the 
unit of energy in our calculations. We have solved the 
system of equations in Eqs. (|33|) and (|34|) for the case of 
the electron hopping-asymmetry in different layers of the 
system. Namely, we suppose the negative value for the 
lower layer -1 hopping parameter £ 1 , associating by this 
the layer-1 with a layer of hole particles. The interlayer 
hopping parameter t± is chosen to be small in order to 
not suppress significantly the exciton gap formation. It 
is clear in Fig. 01 that the excitonic gap is present in the 
case of the absence of the applied external field to the 
system. Thus, A is a spontaneous local gap formed be¬ 
tween the layers in the system and we see in Fig. 0] that 
it is non zero, practically for all values of the in-plane 
Coulomb interaction parameter U/\ti\. The existence of 
the finite gap parameter, at U/\ti\ = 0 , is the fact of 
the dominant role of the interlayer Coulomb interaction 



FIG. 2: (Color online) The interacting band structure, and 
charge gap formation in the bilayer system (see the Ea. (l36l) ). 
The inset shows the symmetry points on a 2D square lattice. 


W. For the intermediate values of the intralayer inter¬ 
action U/\ti\, the gap function endeavor to its maximal 
values. We observe also a specific gap-“see” formation at 
the lower parts of the plots. Indeed, the upper solution 
lines in Fig. 0] represent the strong binding Frenkel-type 
excitons, while the lowest energy solutions, in the form 
of a gap- “see”, represent the all possible Wannier-Mott- 
type excitons in the system. Thus, both type of excitons 
are possible in our bilayer structure in the case of the zero 
applied field to the structure. In Fig.0]the gap functions 
are plotted for the higher value of the interlayer Coulomb 
interaction parameter W/\ti | = 0.75 and also for the 
higher inter-plane hopping amplitude t± = 0.02|£i|. We 
see in Figs. 0] and [5j that there is a smooth passage 
from Frenkel-gap to the Wannier-Mott limit (for exam¬ 
ple at the value U/\ti \ ~ 5.5, for T/|£i| = 0.0001 and for 
w/\h\ ~ 0.75). The merging points “M”, indicated in 
the circles in Figs. 0] and [5j correspond to the points, 
where the two type of solutions are joined. 

As we will see later on, in this Section, the upper so¬ 
lutions for the gap in Figs. 0] and [5] correspond to the 
relatively small values of the chemical potential, while 
the Wannier-Mott gap region is related to the high val¬ 
ues of the chemical potential. In the small and medium 
intralayer interaction regions, both type of excitons coex¬ 
ist. Contrary, for strong intralayer Coulomb interaction 
values, the Frenkel-gap disappears and the fluctuating 
Mott-gap stills open until very high values of the inter¬ 
action parameter U/\ti\. 
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merging point “M”, where the two gap structures join to¬ 
gether. It is clear in Figs. [6] and d that the excitonic gap 
function also increases in the amplitude when increasing 
the interlayer coupling parameter W/|£i|. This fact is 
related to the dominating role of the interlayer coupling 
between the particles. Meanwhile, the parameter t± is 
not affecting considerably the screening effects causing 
by the intralayer Coulomb interaction. 



FIG. 4: Excitonic gap parameter normalized to the elec¬ 
tron layer hopping amplitude, as a function of the in-plane 
Coulomb interaction parameter U/\ti\. Different values of 
the temperature are considered. 


FIG. 3: (Color online) The Monkhorst-Pack grids for the in¬ 
teracting band structure along the high symmetry directions: 
V ^ X, X ^ M, M —T (from left to right). The upper 
panels the grid structure is shown with connected points, and 
in the lower panels the same structure is shown, with uncon¬ 
nected grid-points. 


In Figs. [6] and [7] we have presented the excitonic gap 
parameter evolution for different values of the interlayer 
Coulomb interaction parameter VF/|ti| and for different 
values of the interlayer electron hopping amplitude t±. 
The temperature in all pictures in Figs. [6] and 0is fixed 
at the value T/\ti\ = 0.5. We see in Figs.[6]and [7] that for 
small values of the interlayer Coulomb interaction W/\t\ |, 
the excitonic gaps are very small. Following the evalua¬ 
tions from left to right in Figs. [6] and [7J we observe that 
the excitonic Frenkel-type gap formation peak is displac¬ 
ing to the region of the higher values of the intralayer 
Coulomb interaction parameter U/\ti \ related to the fact 
that for the high values of W/\t \| the screening effects 
of U/\ti\ are relatively small. The same is true for the 


In Fig. [8] we have shown the ecxitonic gap parame¬ 
ter for the case of the hopping-symmetric bilayer, when 
|ti| = |^ 2 1 = 1. The solutions in Fig. [8] are obtained for 
the case W/|ti| = 0.75 and for the interlayer hopping 
amplitude t± = 0.04|£i|. Different values of the tem¬ 
perature are considered. We see in Fig. [8] that there is 
no Frenkel-type solution region and this is related to the 
high hopping value of the conduction band electrons in 
layer-2. Thus, by augmenting the layer-2 hopping am¬ 
plitude, we suppress totally the regular Frenkel-type ex¬ 
citonic gap, and only the fluctuating Wannier-Mott-type 
region remains. 


B. Chemical potential 

In Fig. [9] we have shown the solution for the chemical 
potential for the case W/|ti| = 0.75, t± = 0.04|£i| and 
for T/\t\\ = 0.0001. We see that the solution is given in 
the form of an asymmetric two-band structure, separated 
by the energy gap. Two bands of solutions have opposite 
sign, namely the upper band is positive and the lower 
band has the negative values. The upper band has the 
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FIG. 5: Excitonic gap parameter normalized to the elec¬ 
tron layer hopping amplitude, as a function of the in-plane 
Coulomb interaction parameter U/\ti\. Different values of 
the temperature are considered. 
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FIG. 6: Excitonic gap parameter for small (left panel) and 
intermediate (right panel) values of the interlayer Coulomb 
interaction parameter W/|£i|. The temperature is fixed at 
T/\ti\ = 0.5. 
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FIG. 7: Excitonic gap parameter for intermediate (left panel) 
and large (right panel) values of the interlayer Coulomb in¬ 
teraction parameter VF/|£i|. The temperature is fixed at 
T/\ti | = 0.5. 
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FIG. 8: (Color online) Excitonic gap parameter for the 
hopping-symmetric bilayer: |ti | = 1121 = 1. Different values 
of the temperature are considered. 


sense of the energies, at which the electron creation is 
possible in the layer-2, while the negative band is related 
to the energies of hole-doping in the system, in the layer- 
1. At the intermediate values of the intralayer Coulomb 
interaction U/\ti\ the upper band solutions shoot down, 
to the region of the energy gap, and form a region lo¬ 
cated in the middle of the energy /i-gap. By analyzing 
the solution of /i at T/\ti\ = 0.0001, and corresponding 
values of the excitonic gap parameter A (see in Fig. in, 
we conclude that the solution of fi in the middle of the 
fi- gap is directly related to the degenerated Frenkel-type 
gap in the bilayer structure. 

In Fig.[10]we have presented the solutions for the chem¬ 
ical potential in the large energy scale, which straddle 
all solutions for the case of parameters FF/|£i| = 0.75, 


t l = 0.04|£i| and T/\ti\ = 0.0001. We see that the 
whole solution for the chemical potential forms a leaf¬ 
like structure. This is analogue to the similar structure 
for the chemical potential in the three-dimensional (3D) 
bulk semiconductors and which is discussed in Refsl29l 
and[H. Mainly, the chemical potential, forms a band, in 
3D and 2D semiconducting materials, with the exciton 
formation possibility^ and the upper bound of that so¬ 
lution is responsible for the excitonic gap parameter, and 
for the condensation of the preformed pairs. 

Current discussion concerning the solution for the 
chemical potential could be represented in the schematic 
form given in Fig. [TlJ We see in Fig. [TO] the regions, 
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U/l/,1 


FIG. 9: (Color online) Chemical potential, at T/\ti\ = 0.0001, 
normalized to the hopping amplitude |tij f as a function of the 
in-plane Coulomb interaction parameter U/\ti\. The inter¬ 
layer Coulomb interaction is set as W/\ti \ = 0.75, and the 
interlayer hopping amplitude is t± = 0.04|£i|. 


FIG. 10: (Color online) Chemical potential normalized to the 
hopping amplitude |ti|, as a function of the in-plane Coulomb 
interaction parameter U/\ti\. The interlayer Coulomb inter¬ 
action is set as VF/|£i| = 0.75, and the interlayer hopping 
amplitude is t± = 0.04|£i|. 


where, the Frenkel-type (the red region in the picture) 
or, the Wannier-Mott-type excitons (the region in gray) 
could be formed. The red regions in Fig. [lTJ correspond 
the values of chemical potential responsible for the for¬ 
mation of the Frenkel-type excitons. The red threads 
correspond to the upper solutions of the gap given in 
Figs, [f] and O 

At the merging point “M” (see in Figs. [9] and ITT]) 
the Frenkel-type excitons are “melting” into the see of 
the degenerated Wannier-Mott low-energy region (see the 
low-energy regions in Figs. [4] and [5j). For the strong in- 
tralayer Coulomb interactions U/\ti\, the Frenkel-type 
excitations disappear, and the Wannier-Moot excitons 
are still present and this is corresponding to the van¬ 
ishing of the Frenkel-type gap at the strong values of 
U/\ti\ (see in Figs. H] and The solution region at 

the middle of the fi- gap, which is present in the case of 
very low temperatures (see in Fig. [9] and also the left 
picture in Fig. ID} disappear, when augmenting the tem¬ 
perature (see the right picture in Fig. [TT] and also the 
solutions given in Fig. fT2l) . In Fig. [12] we give the tem¬ 
perature dependence of the chemical potential for the 
case FF/|£i| = 0.75 and t± = 0.04|£i|. Different colors in 
the picture correspond to different values of renormalized 
temperature T/\ti\. The blue points represent the solu¬ 
tions for ii at T/\ti\ = 0.5, red points - for T/\t\\ = 1.0, 
and green points - for T/\ti\ = 1.5. It is clear that the 
temperature kills the solution region in the middle of the 
/i-gap, which is present in the case of very low tempera¬ 
tures T/\ti \ = 0.0001. We see in In Fig.[l2J that with in- 



> 

u 


FIG. 11: (Color online) Schematic representation for the so¬ 
lution of the chemical potential in the quantum bilayer struc¬ 
ture. The red color represents the solution of the chemical 
potential responsible for the formation of the Frenkel-type 
excitons. The left panel corresponds to the case of zero tem¬ 
perature and the right panel corresponds to higher tempera¬ 
tures. 


creasing the temperature the /i-band separation increases 
and the Frenkel-type solution with the “M”-point is dis¬ 
placed to the region of small intralayer interaction region, 
in accordance with the similar behavior of the merging 
point “M” in the excitonic gap structure in Figs. [4] and 

In Fig. [13] we have presented the temperature depen¬ 
dence of the excitonic gap parameter (both, the Frenkel- 
type, and Wannier-Moot-type solutions are given) for the 
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FIG. 12: (Color online) The temperature dependence of the 
chemical potential in the bilayer structure. The interlayer 
Coulomb interaction parameter is fixed at W/\ti | =0.75 and 
t± = 0.04|ti |. The beginning of the Frenkel-type solution (see 
the “M”-point) is displaced to the smaller intralayer interac¬ 
tion region. 
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FIG. 13: (Color online) The temperature dependence of the 
excitonic gap parameter. Different values of the intralayer 
Coulomb interaction parameter U/\ti \ are considered. 


FIG. 14: (Color online) The enlarged Wannier-Mott-region of the 
bilayer excitonic structure, for different values temperature. The in¬ 
terlayer Coulomb interaction parameter is fixed at VT/|ti| = 1-5 and 
the interlayer hopping amplitude is set as t± = 0.1|ti|- The figures in 
panels a-d, show the very small Wannier-Mott-type gap-lobs, which are 
opening in the system. 
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case W/|£i| =0.75 and small interlayer hopping param¬ 
eter t l = 0.04|ti|. We see in Fig. [13] that, when in¬ 
creasing the intralayer Coulomb interaction parameter, 
the gap function also increases, then for the medium val¬ 
ues of the Coulomb interaction U/\ti\ = 4, it decreases 
and finally disappears at the large values ofU/\ti\. This 
behavior is consistent with the behavior of the excitonic 
gap function given in Figs. [4] and [5] 

Furthermore, in the Fig. [14] we show how very small 
Wannier-Moot gap lobs are opening in the spectrum of 
the excitonic gap function. For this, we consider the low¬ 
est gap structure region (see the lowest parts in Figs. 0] 
and [5]). On the large-scale pictures in Fig. [14] we see how 
the very small Wannier-Mott-type gap islands are open¬ 
ing in the gap structure, when increasing the intralayer 
Colomb repulsion. The amplitudes of these small gaps 
are decreasing, when increasing U/\ti\ and vanish, when 
the interlayer Coulomb screening is very high. For the 
very low temperatures (see in the panel-a in Fig. H 
we have only one well defined Mott-gap slope, which is 
spliced in many lobes, when increasing the temperature 
(see the panels-b, -c and -d in Fig. fl4l) . 

V. FINAL REMARKS AND CONCLUSIONS 

We have considered the electronic double layer struc¬ 
ture without the applied external electric field. Employ¬ 
ing the effective action formalism we have shown the com¬ 
plicated excitonic gap structure in the bilayer system. We 
demonstrated that both Frenkel and Wannier-Mott-type 
gaps are possible in the strongly correlated bilayer sys¬ 
tem. Considering the chemical potential structure, we 


are able to separate all possible gap formation regions 
and discuss the related energy scales. For the Frenkel- 
type gap, we have found a very narrow region in the 
chemical potential structure, where the Frenkel-type ex- 
citons could form. 

Meanwhile, we have shown that the Frenkel-type exci- 
tons correspond to the very small values of the chemical 
potential, contrary to the case of the Wannier-Mott low 
excitation region, where the formed excitons correspond 
to the very large particle /i-energy scales. 

The results show that the single-particle energies at 
which the Wannier-Mott-type excitons could be formed 
are relatively very high, and, energetically, it appears 
that the Frenkel-type solution is more favored for the 
bilayer structure. The calculated [/-dependence of the 
excitonic gap function shows that for the small and 
medium values of the Coulomb interaction parameter 
U/\ti \ both, Frenkel and Wannier-Mott-type exciton for¬ 
mations are possible, while in the strong intralayer in¬ 
teraction case the Frenkel-gap is vanishing rapidly, while 
the Wannier-Mott gap survives at relatively high values 
of the Coulomb interaction parameter U/\ti\. For the fu¬ 
ture, it would be interesting to consider the double lay¬ 
ered structure in the case of the non-zero external elec¬ 
tric field and to see the modifications that occur in the 
gap structure for that case. We suppose that the non¬ 
zero applied electric field will favor the exciton formation 
processes in the bilayer system and we expect that the 
amplitudes of the excitonic gap parameters should be 
much higher in that case. Furthermore, we are planning 
to discuss the possibility of the exciton condensation in 
the bilayer type systems. 


* Electronic address: v.apinyan@int.pan.wroc.pl 
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